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ABSTRACT 



Data from the Five College Radio Astronomy Observatory CO Mapping Sur- 
vey of the Taurus molecular cloud are combined with extinction data for a sam- 
ple of 292 background field stars to investigate the uptake of CO from the gas 
to icy grain mantles on dust within the cloud. On the assumption that the 
reservoir of CO in the ices is well represented by the combined abundances of 
solid CO and solid CO2 (which forms by oxidation of CO on the dust), we find 
that the total column density (gas + solid) correlates tightly with visual ex- 
tinction (Ay) over the range 5 < Ay < 30 mag, i.e., up to the highest extinc- 
tions covered by our sample. The mean depletion of gas-phase CO, expressed 
as 5 (CO) = iV(CO)i C c/iV(CO)totai, increases monotonically from negligible levels 
for Ay < 5 to ~0.3 at Ay = 10 and ~0.6 at Ay = 30. As these results re- 
fer to line-of-sight averages, they must be considered lower limits to the actual 
depletion at loci deep within the cloud, which may approach unity. We show 
that it is plausible for such high levels of depletion to be reached in dense cores 
on timescales ~ 0.6 Myr, comparable with their expected lifetimes. Dispersal of 
cores during star formation may be effective in maintaining observable levels of 
gaseous CO on the longer timescales estimated for the age of the cloud. 

Subject headings: dust, extinction — ISM: abundances — ISM: individual (Tau- 
rus Dark Cloud) — ISM: molecules 
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1. Introduction 

CO is a vital molecule in interstellar astrophysics for several reasons. For example, CO 
has long been used as a proxy for H2, allowing the distribution of molecular gas to be mapped 
by means of the readily-observable gas-phase emission lines at millimeter wavelengths. As 
a repository for substantial fractions of the available elemental carbon and oxygen, CO also 
plays an important role in the chemical evolution of the interstellar medium (ISM) as a 
reactant in processes that may lead to synthesis of complex molecules. An important issue 
that affects both of these areas of research is the degree to which CO condenses out of the 
gas onto dust grains in molecular clouds. From a theoretical standpoint, it is expected that 
CO will stick efficiently to grain surfaces at the low temperatures (T < 17 K) prevailing in 
dense regions remote from embedded stars, suggesting that virtually all of the CO could be 
removed from the gas on timescales shorter than cloud lifetimes (Leger 1983; Walmsley et al. 
2004) unless there is an efficient desorption mechanism (Roberts et al. 2007; Oberg et al. 
2007). The presence of solid CO has indeed been inferred from density-correlated depletions 
in the observed gas-phase abundance toward cloud cores (e.g., Bacmann et al. 2002), and 
detected directly in some lines of sight by means of its 4.67 /xm absorption feature in the 
infrared spectra of stars observed through molecular-cloud material (e.g., Chiar et al. 1995; 
Pontoppidan et al. 2003). 

To quantify the distribution of CO between gaseous and solid phases reliably as a func- 
tion of physical conditions is thus an important astrophysical goal. It may be accomplished 
by measuring and comparing gas-phase and solid-state CO column densities in the same 
lines of sight. However, this apparently simple task proves to be quite difficult to accomplish 
in practice. CO-bearing ice mantles that form on the dust are detectable only in absorption, 
by means of their vibrational transitions in the infrared, so observations of the solid phase 
are restricted to lines of sight toward background continuum sources, typically young stellar 
objects (YSOs) embedded within molecular clouds, or field stars located behind them. In 
principle, gas-phase CO may be studied in either absorption or emission, via its electronic 
transitions, vibration- rotation bands or pure rotation features at ultraviolet, infrared and 
millimeter wavelengths, respectively. The most reliable method would be to study solid 
and gaseous phases in absorption simultaneously, in the same lines of sight, against the same 
continuum sources. Acquisition of the necessary data is relatively straight-forward for COi ce , 
far more challenging for CO gas . Studies of the ultraviolet absorption features of gaseous CO 
(e.g., Federman et al. 1994, 2003) are limited to lines of sight with extinctions (Ay < 2) well 
below the level at which ice is detected. The gas-phase lines within the infrared vibration- 
rotation bands are narrow and suffer from significant telluric contamination in ground-based 
observations, requiring both high spectral resolving power (R > 10 4 ) and excellent signal- 
to-noise for effective observation. Data of sufficient quality have been acquired to date only 
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for a small number of relatively bright sources (Mitchell et al. 1988, 1990; Shuping et al. 
2001; Rettig et al. 2005), and these turn out to be exclusively YSOs of high or intermediate 
mass. The further possibility to observe pure rotational transitions in absorption proves 
to be infeasible because of the lack of sufficiently strong background sources at millimeter 
wavelengths and the need to observe multiple transitions to account for the total column 
density, which is spread over many J-states. 

In summary, to date, no convincing detection of interstellar gaseous CO has has been 
made in absorption toward any field star located behind a significant column of dense 
molecular-cloud material, which is unfortunate, as it is the field stars that yield the most 
reliable information on CO depletion in quiescent regions of the clouds, whereas embedded 
YSOs may drive sublimation of ices in their local environment. To measure CO depletion in 
quiescent clouds it is necessary to estimate gas phase column densities from millimeter-wave 
observations of CO emission in lines of sight toward background field stars, for comparison 
with the infrared absorption-line data for solid CO. This raises an obvious concern that the 
two types of observation might not sample identical regions of space, because of differences 
in effective beam size and the possibility of material behind the star, which contributes to 
the millimeter-wave intensity but not to the ice absorption. The well-known problem of 
saturation in the gas-phase CO emission lines is a further concern, generally necessitating 
observation of rarer isotopologues such as 13 CO, C 18 0, C 17 and 13 C 18 (e.g., Frerking et al. 
1982; Bensch et al. 2001a; Harjunpaa et al. 2004). 

This paper compares gas-phase CO column densities extracted from the Five College 
Radio Astronomy Observatory (FCRAO) CO Mapping Survey of the Taurus molecular cloud 
(Goldsmith et al. 2008; Narayanan et al. 2008; see Section 2) with extinction data for some 
292 background field stars in the same Galactic region (Shenoy et al. 2008; Whittet et al. 
2001). The spatial resolution of the gas-phase CO observations from the FCRAO survey 
(~ 45") is improved by a factor ~ 2 compared with an earlier study by Frerking et al. 
(1982) for a sample of 14 background stars in the Taurus region. Contamination of the 
gas-phase column densities by molecular material behind the stellar sources in our sample is 
unlikely to be a significant problem, given the location of the cloud away from the Galactic 
equator (b ~ 15°) in a direction where virtually all of the observed interstellar material is 
confined to the cloud itself. Column densities of the ices have been shown in previous work 
on a subset of our sample to be well correlated with Ay (Whittet et al. 2007 and references 
therein), enabling reasonable estimates to be made in cases where direct measurements are 
unavailable. Gas-phase and solid-state abundances may thus be compared over the entire 
sample. We use our results to show (Section 3) that the mean line-of-sight depletion of 
CO from gas onto dust increases systematically from negligible levels at low extinction to 
~ 60% at Ay ~ 30. Our results are in broad agreement with theoretical predictions based 
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on expected timescales for the depletion of gaseous CO in quiescent cloud cores (Section 4). 



2. Gas-phase CO column density determination 

We use the 12 CO and 13 CO observations of the Taurus molecular cloud presented by 
Goldsmith et al. (2008) to obtain gas-phase column densities, iV(CO) gas , at the positions of 
our program stars. The primary sample is a catalog of 247 reddened background field stars 
identified from Two-Micron All Sky Survey (2MASS) photometry and other infrared data 
by Shenoy et al. (2008, their Table 1), covering the extinction range 3 < Ay < 29. To enable 
better coverage of low extinctions, we add optically-selected field stars with well-determined 
extinction values from Whittet et al. (2001) and references therein (see Figure 5 of Shenoy 
et al. 2008 for a histogram of the distributions), thus increasing the sample size to 292. 
Finally, we include for comparison ten highly reddened YSOs from the Shenoy et al. study 
(their Table 2). All stellar positions used in the current work are taken from the 2MASS 
Point Source Catalog (Skrutskie et al. 2006). The complete sample is listed in Table 1 (online 
version) . 

We first calculated the column density of 13 CO assuming that its J = 1 — > transition 
is optically thin, that 12 CO J = 1 — > is optically thick, and that local thermodynamic 
equilibrium applies. The column density of 13 CO is then proportional to its integrated 
intensity and inversely to the spontaneous decay rate. It also depends on the excitation 
temperature, derived from 12 CO, which affects the upper level population and the partition 
function. We also determine the optical depth, which is a function of the ratio between the 
13 CO intensity and the excitation temperature, to apply a correction for saturation. Finally, 
we transform the column density of 13 CO to that of 12 CO assuming an isotopic ratio of 65 
(Bensch et al. 2001a) in the relatively well-shielded portions of the cloud defined by the 
"mask 2" region of Goldsmith et al. (2008). 

The determination of A^(CO) gas was improved relative to that presented by Goldsmith 
et al. (2008) by including an updated value of the spontaneous decay rate, using an exact 
numerical rather than approximate analytical calculation of the partition function. Ad- 
ditionally, the data were corrected for error beam pick-up using the method presented by 
Bensch et al. (2001b). The resulting values of iV(CO) gas are about ~20% larger than those of 
Goldsmith et al. (2008). Full details of these improvements are presented elsewhere (Pineda 
et al. 2010). 

Motivated by observations of core-to-edge temperature differences in molecular clouds 
(Evans et al. 2001) which can be found even in regions of only moderate radiation field inten- 
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sity, Pineda et al. (2010) studied the effects of temperature gradients on the determination 
of iV(CO) gas . The radiative transfer code RATRAN (Hogerheijde & van der Tak 2000) was 
used to model the 12 CO and 13 CO emission emerging from a model cloud. Pineda et al. 
found that using 12 CO to determine the excitation temperature of the CO gas only traces 
the temperature at low column densities while the excitation temperature is overestimated 
for larger column densities. This produces an underestimate of the 13 CO optical depth, 
and in consequence the opacity correction of iV( 13 CO), which is usually evaluated assuming 
an isothermal cloud. The column densities presented here were corrected by this method 
including modest edge-center temperature gradients of ~ 4 K. This procedure typically in- 
creases the value of iV(CO) gas by ~10-40%, with larger increases occurring at higher column 
densities. Final values are presented in Table 1, together with ice and extinction data. 



3. Results and discussion 

3.1. Gaseous CO versus extinction 

A plot of N(CO) gas against Ay is shown in Figure 1. In addition to results from the 
current work, we include A^(CO)v as estimates from the data of Frerking et al. (1982) for lines 
of sight to 14 Taurus field stard_| in combination with Ay values from Shenoy et al. (2008). 
Previous observations of gas-phase C 18 toward several molecular clouds suggest a general 
linear relation 

iV(C 18 0) gas = 2 x 10 14 (A^ - 2) cm" 2 (1) 

for Ay values up to about 10 mag (see Kainulainen et al. 2006 and references therein), 
and this relation is represented (after conversion to the common isotopic form) by a dotted 
line in Figure 1. The distribution of the data is, indeed, consistent with this relation at 
low extinction. However, a clear divergence is apparent for Ay > 10, with little further 
increase in iV(CO) gas with Ay. The general trend is well represented by a sigmoidal fito 
shown as a solid curve in Figure 1. It is possible that this behavior might result in part from 
residual saturation in our column density calculations (Section 2), leading to underestimates 
at the highest optical depths; however, similar behavior is seen in the data of Frerking et al. 
(1982) for rarer isotopologues such as C 17 and 13 C 18 that are presumed to be unsaturated 



lr The abundance ratios we used to estimate -ZV(CO) gas in the common isotopic form from the various 
isotopologues observed by Frerking et al. (1982) were 12 CO/ 13 CO = 65, C 16 0/C 18 = 560, and C 16 0/C 17 
= 1790 (Bensch et al. 2001a; Ladd 2004; Wilson & Rood 1994). 

2 The functional form is y = a 2 + (o&i — a,2)j{\ + (x/xq) p ); see |http://www.originlab.com Note that the 



fit is purely empirical, for illustrative purposes, and has no basis in theory. 
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(Figure 1), and we therefore consider residual saturation in our data to be at worst a second 
order effect. Note also that a similar trend in iV(CO) gas with Ay was found in another 
quiescent dense cloud (IC 5146) by Kramer et al. (1999), based on observations of C 18 0. 



3.2. Adding the solid phase 

Previous observations have shown that both CO and CO2 are present in the ices over 
the same range of Ay in which gas-phase CO appears to become depleted in the Taurus 
cloud (Whittet et al. 1989, 1998, 2007; Chiar et al. 1995; Bergin et al. 2005; Knez et al. 
2005). Astrochemical models and laboratory simulations predict that CO2 forms in situ on 
interstellar grains, as the result of simple surface oxidation reactions such as CO + O — > CO2 
and CO + OH ->■ C0 2 + H (Tielens & Hagen 1982; Ruffle & Herbst 2001; Roser et al. 2001). 
A substantial fraction of the CO that becomes depleted onto dust is thus predicted to be 
converted to C0 2 on timescales short compared with cloud lifetimes, and so a complete 
inventory of adsorbed CO should include the reservoir represented by solid CO2 as well as 
by CO itselo An estimate of the "total CO" column density that accounts for both depleted 
CO and that remaining in the gas is therefore given by 

iV(CO) total = iV(CO) gas + iV(CO) icc + iV(C0 2 ) ice . (2) 

No other CO-bearing molecule appears to be sufficiently abundant in quiescent clouds to 
constitute a significant reservoir of CO in the solid phase: methanol (CH 3 OH), for example, 
is sometimes abundant toward YSOs but appears to be no more than a minor constituent 
of the ices observed toward field stars (Chiar et al. 1996; Gibb et al. 2004). 

The observations also show that each of the major ice constituents (H 2 0, CO and C0 2 ) 
displays a tight linear correlation between column density and extinction. These correlations 
may be expressed in the form 

N icc = q(A v - A y °) (3) 

where Ay is the "threshold" extinction for detection of a given species, such that A^; ce = for 
Ay < Ay. The best available values of the constants q and Ay for the two species relevant 
to the present study are: 

CO : q = (0.400 ± 0.060) x 10 17 cm" 2 mag -1 , A y ° = 6.7 ± 1.6 mag (4) 

C0 2 : q = (0.252 ± 0.036) x 10 17 cm" 2 mag -1 , A y ° = 4.3 ± 1.0 mag. (5) 



3 Note that CO2 has negligible abundance in the gas-phase in cold molecular clouds: interstellar gaseous 
CO2 is observed only in lines of sight to luminous YSOs that drive rapid sublimation of the ices (van Dishoeck 
et al. 1996). 
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These values are based on fits to data for a sample of 13 Taurus field stars in the extinction 
range 5 < Ay < 24 mag (see, e.g, Figure 4 of Whittet et al. 2007). The majority of stars in 
our sample lack the spectral data that would enable direct measurement of ice column den- 
sities: in these cases, we use Equation (3) with the above values of q and Ay to estimate the 
contributions of solid CO and CO2 to the total column density in Equation (2) from the Ay 
values for field stars. Resulting values are listed in Table 1. We consider the well-established 
iVi ce versus Ay correlations to be sufficiently strong to justify this approach for field stars, 
but it may not be appropriate for YSOs. For this reason, total CO column densities are cal- 
culated for YSOs only in cases where direct measurements of both iV(CO)i ce and iV(C02)i ce 
are available from the literature: only two out of ten in our sample (J04395574+2545020 
and J04400800+2605253) satisfy this requirement (Cook et al. 2010). 

The resulting plot of iV(CO) tota i versus Ay is shown in Figure 2. The overall distribution 
of points is clearly consistent with the linear least-squares fit (solid line), yielding a correlation 
only marginally different from that predicted by Equation (1). This corresponds to an 
abundance Nco/N-r ~ 5.3 x 10~ 5 at extinctions above the threshold, assuming the canonical 
value of Nft/Ay ?» 1.9 x 10 21 cm -2 mag -1 for the total hydrogen gas to extinction ratio 
(Bohlin et al. 1978). 



3.3. The CO depletion factor 

The CO depletion factor is defined 

S(CO) = iV(CO) icc /iV(CO) total = 1 - iV(CO) gas /iV(CO) total , (6) 

such that < 5 (CO) < 1, the lower and upper bounds corresponding to all CO in gaseous 
and solid forms, respectively. This quantity was calculated for each line of sight in our 
sample from the data used to construct Figures 1 and 2, and the results plotted against 
Ay in Figure 3. The curve is calculated in the same way by ratioing the sigmoidal fit from 
Figure 1 with the linear fit from Figure 2. In general, the data in Figure 3 are consistent 
with a monotonic increase in 5 (CO) from zero at low extinction (Ay < 5) to ~ 0.6 at the 
high end of the observed range. The scatter appears somewhat asymmetric because it arises 
primarily in the denominator of Equation (6); outliers well above the trend arise because a 
few lines of sight with intermediate extinction have unexpectedly low gaseous CO column 
densities for their Ay values. 

It is important to bear in mind that 5 (CO) measures the mean depletion along the line 
of sight, and this may differ from the local value at a given location corresponding to a given 
optical depth. In particular, the actual depletion at high Ay deep inside the cloud is expected 



to exceed the value of ~ 0.6 suggested by Figure 3, because of the "dilution" effect of lower- 
density regions that are inevitably included in any line of sight passing through the cloud. 
To illustrate this, we compare data in Table 2 for typical lines of sight at Ay = 5, 10 and 20 
(average values) and at Ay = 28.7 (for the individual star J04135352+2813056, which has 
the highest extinction in our sample). The final column lists differential depletion factors, 
calculated by replacing total line-of-sight column densities in Equation (6) with incremental 
values relative to the previous row: for example, the value 5 (CO) = 0.69 represents a segment 
of column through the cloud extending from Ay = 10 to Ay = 20. Comparing mean and 
differential values of 5 (CO) in Table 2, we find that the latter increase more rapidly with 
Ay, and approach 0.9 (i.e., 90% depletion from the gas) at the highest extinctions sampled 
by our data. 

Our results show that CO depletion is not only high within dense cores but is also 
substantial at intermediate levels of extinction. This is contrary to what is often assumed. 
For example, Padoan et al. (2006) report significant differences in the distributions of gas and 
dust in the Taurus cloud based on a comparison of maps of gas-phase CO and extinction, 
which they attribute to spatial fluctuations in the gas-to-dust ratio driven by turbulent 
flows. These authors assume in their analysis that the CO remains entirely in the gas for 
Ay < 10 mag; however, our results indicate that the depletion is already significant (~ 30% 
or more) at this level of extinction. 



3.4. Comparison of field stars and YSOs 

The ten young stellar objects in our sample display an unremarkable distribution in 
Figure 1, having gaseous CO column densities similar to those of field stars in the same 
extinction range. For reasons discussed above, only two YSOs are included in the calculation 
of total CO column density and depletion (Figures 2 and 3). These objects have estimated 
mean depletion factors of 0.36 (J04400800+2605253) and 0.47 (J04395574+2545020), values 
that are again consistent with those of field stars in the same extinction range. This is 
in contrast to the results of Mitchell et al. (1990), who find 5 (CO) < 0.1 for a sample of 
seven YSOs. The most likely explanation is a systematic difference in mass between the two 
samples: massive YSOs such as those in the Mitchell et al. study evolve rapidly and become 
sufficiently luminous to drive widespread sublimation of interstellar ices in their vicinity, 
whereas the lower-mass objects in our sample have less impact on their environment. For a 
young star of a given mass, the degree of depletion of CO and other molecules in the envelope 
is likely to evolve with time (e.g., Thomas & Fuller 2008). Rettig et al. (2005) estimate that 
the solid CO abundance exceeds the gas phase abundance by a factor ~ 6 toward the eruptive 
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T Tauri star V1647 Orionis, which implies 5 (CO) > 0.85 in this line of sight. It would clearly 
be important to obtain data for a larger sample of low and intermediate-mass YSOs. 

In the case of one object in our sample, the intermediate- mass young star Elias 18 
(J04395574+2545020), a notable discrepancy is apparent between our value of iV(CO) gas ~ 
1.0 x 10 18 cm~ 2 from the radio survey data and that of Shuping et al. (2001), who obtain a 
value a factor ~ 100 less from a study of the infrared spectrum of this source. It is possible 
that the difference might arise because of the finite beam size of the radio observations, and 
the possibility of dense material behind the source, as previously discussed in Section 1. 
Shuping et al. propose that CO is highly depleted in the circumstellar disk of Elias 18 and 
that it suffers very little foreground absorption in the molecular cloud; however, we consider 
this unlikely, given that the abundance of solid CO is actually less toward Elias 18 compared 
with the Taurus field stars with respect to both H 2 0-ice and total extinction (Ay « 22 mag) 
in the line of sight (Chiar et al. 1995; Nummelin et al. 2001). Hence, the total CO as well 
as the gas-phase CO toward Elias 18 must be unusually low if the Shupping et al. result is 
accurate. 



4. Comparison with models 

In this section we briefly review and discuss theoretical constraints on the timescale 
for CO depletion in molecular clouds, and examine the degree to which they are consistent 
with our results. The rate at which gaseous atoms or molecules accumulate onto the dust 
depends on the collision rate and the probability that a collision will lead to attachment. 
The collision rate depends on the thermal speed of the gaseous species and the surface area 
of the dust. The latter may be expressed in terms of the total grain geometric cross- sectional 
area per H atom: 

Ed = f— r„(q)Mo.) da (7) 

/„ . % 

•J <2 mm n 

where a represents an appropriate grain dimension such as the radius for a sphere, o~&(a) is the 
cross-sectional area of a grain, n^(a) is the number density of dust in the size range a — > a+da, 
and nu = n(HI) + 2n(B.2) is the total hydrogen (proton) density (see, e.g., Spitzer 1978; 
Draine 1985). It is usual to evaluate S^ assuming spherical particles and a power-law size 
distribution n^(a) oc a -3 ' 5 (the "MRN" size distribution of Mathis et al. 1977), with the upper 
and lower limits a max and a min set by optimization of grain models fit to the mean interstellar 
extinction curve for the diffuse ISM. This approach yields Ed ~ 1.0 x 10~ 21 cm 2 H _1 . 

Assuming that grain speeds are small compared with the mean thermal speed (v) of 
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impinging gas particles, the timescale for depletion is given by 

^dep = rr; — (8) 

where v = (SkT^/jcm) 1 ' 2 for particles of mass m, and £ is the sticking probability. Both 

theoretical models and experimental data predict a high sticking probability (£ ~ 1) for CO 

at molecular-cloud temperatures (Leger 1983; Bisschop et al. 2006). Assuming T gas = 12 K 

(e.g., Stepnik et al. 2003) and setting m to the mass of CO, Equation (8) then yields the 

result 

3 x 10 9 

tdep ~ yr (9) 

n H 

where % is in units of cm -3 . 

A systematic error in estimates of tdcp that is generally overlooked arises because the 
calculation of S d assumes a grain size distribution appropriate to diffuse (low-density) phases 
of the ISM. In dense clouds, the grains grow by coagulation as well as by adsorption of 
atoms and molecules from the gas. This can lead to a substantial reduction in the relative 
number of small grains, and hence in cross-sectional area, on timescales short compared with 
cloud lifetimes. Draine (1985) estimates that for spherical grains following the MRN size 
distribution, 80% of Ed is contributed by grains with radii a < 0.05 /zm that contain only 
~30% of the total mass of dust. Coagulation may thus lead to substantial reduction in 
Ed, with proportionate increase in tdcp- Observational evidence for growth by coagulation 
is provided by systematic changes in the extinction curve, characterized by the parameter 
Ry = Av/Eb-v (the ratio of total to selective extinction; Fitzpatrick 1999; Whittet 2003). 
Consider, for example, ultraviolet extinction at wavelength A = 0.1 /im, which is where 
particles of radius a ~ A/2-7T ~ 0.02 /im contribute most efficiently: the extinction is observed 
to decline by a factor ~2 as Ry increases from the average diffuse-ISM value of 3.1 to a 
typical value of 5.5 within a dense cloud (Whittet et al. 2001), implying a corresponding 
reduction in n^ for grains in this size range. Equation (9) may thus underestimate the 
depletion timescale within molecular clouds by a factor of order 2. A typical dense core 
within the Taurus complex has density nn ~ 10 4 cm -3 (see, e.g., Table 1 of Onishi et al. 
1996). Allowing for a doubling of the constant in Equation (9), this density corresponds to 
a CO depletion timescale ~ 0.6 Myr. 

The age of the Taurus complex is at least a few Myr (e.g., Kenyon & Hartmann 1995). 
Estimates of this order for Taurus and other dark clouds are often used in the literature 
in conjunction with variants on Equation (9) to support the view that all the available CO 
is expected to be in solid form (and hence to be unobservable at millimeter wavelengths) 
on timescales much less than cloud lifetimes, unless an efficient (and poorly understood) 
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desorption mechanism is operative (e.g., Roberts et al. 2007, Oberg et al. 2007 and references 
therein). For a cloud such as Taurus, in which star formation has been ongoing for at least 
1-2 Myr (Kenyon & Hartmann 1995), the mechanism may simply be star formation itself. 
Protostars form within prestellar cores and disperse them on timescales < 1 Myr (Kenyon & 
Hartmann 1995; Onishi et al. 2002; Visser et al. 2002). The life expectancy of an individual 
core is a factor of 5-10 less than that of the entire complex. 

To summarize, in general these estimates seem consistent with our results. It is highly 
plausible that CO depletion factors exceeding 50% can be reached within dense cores on 
timescales comparable with their expected lifetimes. However, no crisis is implied in our 
understanding of desorption mechanisms by the observation that gas phase CO remains 
observable over the range of densities covered by our survey. 



5. Conclusion 

We have studied the uptake of CO from gas to dust in the Taurus molecular cloud by 
comparing CO column densities from millimeter-wave observations with extinction values 
and ice column densities determined from infrared observations of stars. The main conclu- 
sion of this work is that the mean depletion of gas-phase CO increases monotonically from 
negligible levels at visual extinctions Ay < 5 to about 60% at Ay ~ 30. The depletion is 
significant (~30%) even at moderate levels of extinction (Ay ~ 10). These results refer 
to line-of-sight averages and should be considered as lower limits to the actual depletion 
at loci deep within the cloud, which may approach 90-100% in dense cores. We estimate 
the timescale for CO depletion in a core of density % ~ 10 4 cm -3 to be ~ 0.6 Myr, a 
value consistent with independent estimates of the lifetime of a typical core. Dispersal of 
cores during star formation may help to maintain observable levels of gaseous CO over cloud 
lifetimes. Our results will provide useful constraints on astrochemical models for molecular 
clouds (e.g., Chang et al. 2007; Cuppen & Herbst 2007), and on estimates of macroscopic 
cloud properties such as mass, structure, and the distribution of molecular material (e.g., 
Padoan et al. 2006; Pineda et al. 2010). 

The constraints placed on CO depletion by our analysis probably represent the limits 
of what can be accomplished by combining data from the disparate techniques of gas-phase 
emission spectroscopy and solid-state absorption spectroscopy. Further refinement awaits a 
targeted survey of infrared vibration-rotation lines in absorption toward background field 
stars known to display ice features, requiring 4.5-4.8 /im spectra of high signal-to-noise with 
resolving powers > 10 4 . We also advocate further study of intermediate-mass YSOs by this 
technique, to investigate the degree to which they may contribute to the global recycling 
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of CO between solid and gaseous phases in molecular clouds, and to better evaluate other 
astrochemical pathways for CO, such as CH 3 OH formation, in protostellar envelopes. 
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trobiology Institute (grant NNA09DA80A). This work was supported in part by the Jet 
Propulsion Laboratory, California Institute of Technology, under a contract with NASA. We 
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Table 1. Extinction values and column densities for program stars a . 



n 


Identification 


A v 


7V(CO) gas 


N(CO) icc 


iV(C0 2 )i C0 


[1] 


[2] 


m 


[4] 


[5] 


[6] 


1 


J04082673+2803429 


6.4 


4.91 


0.00 


0.53 


2 


J04090144+2453214 


1.4 


3.14 


0.00 


0.03 


3 


J04092063+2816031 


1.1 


5.10 


0.00 


0.00 


1 


J04104300+2820340 


4.6 


1.99 


0.00 


0.08 


5 


J041 10488+2443185 


3.5 


2.30 


0.00 


0.00 


6 


J041 12677+2831093 


5.5 


7.67 


0.00 


0.30 


7 


J04112801+2830271 


3.7 


6.91 


0.00 


0.00 


8 


J04113168+2829562 


1.5 


6.69 


0.00 


0.05 


9 


J04114185+2841282 


5.9 


2.51 


0.00 


0.40 


10 


J04114938+2815568 


3.6 


4.92 


0.00 


0.00 


11 


J04122296+2949441 


5.6 


2.85 


0.00 


0.33 


12 


J04123318+2945152 


1.8 


4.75 


0.00 


0.13 


13 


J04123940+2816419 


7.0 


6.64 


0.12 


0.68 


11 


J04130664+2235365 


5.2 


1.88 


0.00 


0.23 


15 


J04132688+2804584 


6.3 


3.52 


0.00 


0.50 



a The complete table is available in the online version of this paper. 
Description of columns: [1] Running number. [2] Source identification: field stars from 
Shenoy et al. 2008 (n = 1 — 247) designated by 2MASS catalog number; additional 
field stars (n = 248 — 292) mostly designated by Henry Draper (HD) catalog number 
(Whittct ct al. 2001); YSOs (n = 293-302) designated by 2MASS catalog number. 
[3] Visual extinction in magnitudes. [4] Gas phase CO column density from the present 
work in units of 10 17 cm -2 . [5] and [6] Solid phase CO and CO2 column densities, 
respectively, in units of 10 17 cm -2 ; direct measurements are included where available 
(Whittet et al. 2007, their table 2, for field stars; Cook ct al 2010, their table 3 for 
YSOs), all other values are calculated using the method described in Section 3.2. 
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Table 2. Mean and differential CO depletion factors for field stars. 



A a 


Af(CO) gas 


iV(CO) tota i 


<5(CO) 6 


<5(CO) c 


(mag) 


(xl0 17 cm~ 2 ) 


(xl0 17 cm~ 2 ) 


(mean) 


(difference) 


5 


.3.9 


4.2 


0.07 


— 


10 


7.0 


9.6 


0.27 


0.43 


20 


10.2 


19.9 


0.49 


0.69 


28.7 


10.9 


25.8 


0.58 


0.88 



a Data at Ay = 5, 10 and 20 mag arc averages for typical lines of sight 
falling within 20% of these extinction values; data at Ay = 28.7 mag 
refer to the individual star J04135352+2813056. 

b Values calculated (Equation 6) from line-of-sight column densities. 

c Valucs calculated (Equation 6) from incremental differences in col- 
umn densities (current row relative to previous row). 
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Fig. 1. — Plot of gaseous CO column density against visual extinction for lines of sight in our 
sample (Table 1). Crosses (x) and star symbols represent background field stars and young stellar 
objects selected from the 2MASS Point Source Catalog (Shenoy et al. 2008), respectively. Pluses 
(+) near the origin represent additional optically-selected field stars (Whittet et al. 2001). Also 
plotted are data based on gas-phase 13 C 18 0, C 17 and C 18 observations toward selected field 
stars by Frerking et al. (1982), represented by diamonds, squares and circles, respectively. The 
diagonal dotted line represents the empirical relationship N(C ls O) = 2 x 10 14 (A\/ — 2) cm -2 based 
on data for several molecular clouds (Kainulainen et al. 2006 and references therein), assuming 
N(CO)/N(C ls O) = 560. The solid curve is a sigmoidal fit to the data for the entire sample (see 
Section 3.1; the fit parameters are a\ = —1.532, 02 = 23.202, xq = 22.541 and p = 0.8164). 
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Fig. 2. — Plot of total (gas + ice) CO column density against visual extinction. Solid CO2 is 
also included in the ice component (see text, section 3.2). Symbols have the same meaning as in 
Figure 1. The dotted line represents the empirical relationship iV(C 18 0) = 2 x 10 14 (Ay — 2) cm~ 2 , 
as in Figure 1. The solid line is the unweighted linear least-squares fit to the data for the entire 
sample, given by 7V( C0 )totai = (1-024 ± 0.025) [A v - (0.78 ± 0.15)] x 10 17 cm' 2 . 
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Fig. 3. — Plot of mean CO depletion factor 5 (CO), denned in Equation (6), against visual ex- 
tinction. Symbols have the same meaning as in Figures 1 and 2. The solid curve is obtained by 
combining the sigmoidal fit from Figure 1 with the linear fit from Figure 2. 



